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ABSTRACT: With the large-scale integration of renewable en- 
ergy and grid-following power electronic devices, system volt- 
age support strength decreases, thus threatening the system's safe 
and stable operation. In the homogeneous scenario where power 
electronic devices integrated into the system have the same dy- 
namics, a theoretically rigorous and highly stable quantitative 
method for quantifying system voltage support strength can be 
formed based on the generalized short-circuit ratio(gSCR) and 
the device or station critical short-circuit ratio(SCRo) In the het- 
erogeneous scenario where power electronic devices integrated 
into the system have weakly different dynamics, system voltage 
support strength can be quantified by the first-order approxima- 
tion of gSCR and SCRo, based on some special characteristics of 
devices and power grid. However, there is a lack of unified der- 
ivation principles and calculation methods. To this end, this pa- 
per focuses on quantifying voltage support strength under small- 
signal stability and discovers that the multi-infeed system can be 
approximately decoupled into multiple low-dimensional sys- 
tems. On this basis, this paper proposes the concept of eigen- 
subsystem and its calculation method and interprets their physi- 
cal significance. Based on the concept of eigen-subsystems, the 
derivation principle and calculation method of gSCR and SCRo 
are given. Additionally, the specific calculation methods for in- 
dices are provided when considering grid-following converters 
under non-rated operating conditions, reverse active power out- 
put and grid-forming devices. Finally, the effectiveness of the 
principles and methods is verified in several cases. 
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摘要 : 随 着 跟 网 型 新 能 源 /电力 电子 装备 的 大 规模 接 入 ， 系 
统 电压 支撑 强度 降低 ， 新 型 电力 系统 的 安全 稳定 风险 增加 。 
在 同 构 新 能 源 /电力 电子 装备 接 入 场景 下 ， 利 用 电网 广义 短 
路 比 和 装备 / 场 站 临界 短路 比 可 以 形成 与 稳定 性 强 相 关 的 系 
统 电压 支撑 强度 量化 方法 ; 在 弱 异 构 场 景 下 , 基于 装备 和 电 
网 动态 的 特殊 性 质 ， 利 用 广义 短路 比 和 装备 / 场 站 临界 短路 
比 的 一 阶 近似 可 以 实现 强度 量化 , 但 缺乏 统一 的 导出 原理 及 
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计算 方法 。 为 此 , 本 文 聚焦 小 扰动 维度 下 电压 支撑 强度 量化 
问题 , 首先 利用 多 馈 入 系统 可 近似 解 耦 为 多 个 低 维 子 系统 的 
规律 , 提出 特征 子 系统 的 概念 和 计算 方法 , 并 诠释 了 其 物理 
意义 ; 其 次 , 基于 特征 子 系统 , 提出 了 广义 短路 比 及 装备 / 场 
站 临界 短路 比 的 导出 原理 及 计算 方法 。 此 外 , 针对 电力 电子 
装备 在 非 额 定 运行 点 、 部 分 装备 有 功 功 率 反 向 以 及 包含 构 网 
型 装备 的 场景 下 ,给 出 了 广义 短路 比 的 具体 计算 公式 。 最 后 ， 
算 例 验 证 了 所 提 原 理 和 方法 的 有 效 性 。 


关键 词 : 广义 短路 比 ， 装 备 / 场 站 临界 短路 比 ， 特 征 子 系统 
0 引言 


随 着 新 能 源 的 快速 发 展 ， 电 力 系 统 逐 步 演化 成 
含有 高 比例 新 能 源 及 电力 电子 装备 (后 文 简称 为 
“装备 ”) 的 新 型 电力 系统 中 中。 电压 支撑 强度 表征 
了 扰动 下 系统 电压 抵抗 偏 移 以 及 抵抗 失 稳 的 能 力 ， 
强度 高 表明 系统 抵御 被 外 力 破坏 的 能 力 强 ， 受 扰 后 
电压 响应 性 能 好 。 然 而 ， 具 有 强 支 撑 作 用 的 同步 机 
占 比 减少 ， 新 能 源 集 中 并 网 容量 大 、 与 主 网 电气 距 
离 远 且 常 采用 跟 网 型 控制 策略 而 对 电网 支撑 强度 
弱 ， 使 得 新 型 电力 系统 电压 支撑 强度 变 低 ， 受 扰 后 
安全 稳定 风险 增加 B77， 吸 需 科 学 的 强度 量化 和 计 
算 方法 (考虑 到 低 惯 量 维度 的 强度 问题 可 以 分 开 考 
谍 ， 后 文 为 表述 方便 ,将 系统 电压 支撑 强度 简写 为 
“系统 强度 ”)。 

大 规模 跟 网 型 装备 (如 不 特别 说 明 ， 后 文 简称 
“装备 ”) 接 入 后 , 系统 强度 的 分 析 计 算 存 在 以 下 难 
点 : 装备 及 电网 母线 数量 多 ， 且 装备 与 网 络 动态 交 
互 复杂 ， 高 维 动态 特征 难以 简化 ， 装 备 种 类 多 ， 如 
空 制 策 略 和 参数 等 各 异 的 跟 网 型 变 流 器 或 构 网 型 
变 流 器 等 BTINl， 进 一 步 增加 了 系统 的 复杂 度 ; 新 能 
源 运行 工 况 多 ， 也 给 强度 量化 带 来 巨大 挑战 。 

为 分 析 系 统 强 度 且 方便 工程 应 用 ， 工 业界 常用 
电网 强度 近似 表征 系统 强度 ， 并 引入 单 馈 入 和 多 人 馈 
入 短路 比 作 为 量化 指标 。 其 中 ， 电 网 强度 表征 了 在 
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不 计 及 跟 网 型 装备 动态 作用 下 母线 电压 在 扰动 作 
下 抵抗 偏离 正常 运行 点 的 能 力 。 例如 : CIGRE 多 
馈 入 短路 比 (multi-infeed short-circuit ratio, 
MSCR)H31, 位置 相 关 短 路 比 (site-dependent SCR, 
SDSCR)!!41, 4 Yq wh Fae LL (multi-renewable energy 
station SCR, MRSCRJ)U5、 等 效 短路 比 (equivalent cir- 
cuit-based SCR，ESCRJ)023 和 带 交 互 因子 的 短路 比 
(the SCR with interaction factors, SCRIFP)D] 等 指标 。 
上 述 指标 可 以 反映 装备 接 入 前 后 母线 电压 稳 态 1 
移 量 , 但 由 于 未 考虑 装备 间 及 其 与 电网 的 宽频 耦合 
特性 ， 故 与 系统 稳定 性 关联 弱 , 用 于 反映 稳定 特性 
I oR ZS LER SCTE 

文献 [7]、[16]-[18] 从 装备 和 电网 动态 交互 视角 
E 塑 了 融合 装备 和 电网 动态 的 系统 强度 内 涵 ， 并 提 
基于 电压 对 电流 (或 电压 对 功率 ) 灵敏 度 方程 的 
电网 广义 短路 比 (generalized short-circuit ratio, gSCR) 
指标 和 装备 / 场 站 临界 值 指标 , 从 理论 上 实现 了 系统 
强度 与 小 扰动 稳定 性 的 关联 。 在 同 构 系 统 〈 即 接 入 
装备 的 动态 特性 相似 ， 有 具体 定义 见 文 献 [7]) 中， 多 
电力 电子 装备 馈 入 系统 (简称 “多 馈 入 系统 ”) 可 解 
耦 为 装备 动态 整体 保持 的 多 个 等 效 单 馈 入 系统 ; 基 
等 效 单 馈 入 系统 ， 类 比 传统 短路 比 的 用 途 ， 定 义 
电网 广义 短路 比 和 装备 / 场 站 临界 短路 比 , 并 论证 
了 装备 / 场 站 临界 短路 比 等 于 广义 短路 比 临界 值 的 
寺 论 ， 从 而 实现 了 考虑 装备 动态 特性 的 系统 强度 量 
; 在 弱 异 构 场 景 下 ， 如 外 特性 存在 差异 的 直 驱 风 
机 上 、 光伏 本 、 构 网 型 装备 0000 和 储 能 系统 接 入 09] 
或 者 装备 在 非 额定 运行 点 [等 场景 , 异 构 装备 动态 
可 看 作 是 同 构 系统 中 装备 或 网 络 动 态 的 扰动 ， 利 用 
广义 短路 比 和 装备 / 场 站 临界 短路 比 的 一 阶 近似 即 
可 实现 强度 量化 。 上 述 广 义 短路 比 及 其 临界 值 的 计 
算 关 键 为 求解 能 够 近似 原 系统 关键 模 态 稳定 性 的 
解 耦 单 馈 入 系统 。 然 而 ， 在 这 些 非 理想 场景 中 ， 寻 
找 等 效 解 耦 单 馈 入 系统 的 理论 和 方法 缺乏 普 适 性 ， 
不 利于 广义 短路 比 在 复杂 场景 下 的 应 用 。 

为 此 , 本 文 则 在 统一 多 种 场景 下 广义 短路 比 的 
导出 原理 及 计算 方法 。 首先 , 引入 特征 子 系统 的 概 
念 并 提出 其 求解 方法 ; 其 次 , 基于 特征 子 系统 提出 
广义 短路 比 与 装备 / 场 站 临界 短路 比 的 通用 导出 原 
理 和 计算 方法 ; 最 后 , 算 例 验证 了 方法 的 有 效 性 。 
1 预备 知识 
1.1 FERRET SST CREA 

2nx2n FERE M 的 特征 子 空间 及 广义 特征 矩阵 

定义 如 下 : 

定义 120〈 特 征 子 空间 和 广义 特征 矩阵 ): 记 
R(X) ={Xuu € C” } A R(Y ) = {Yuu €C” } 4H 


=r 


m. cr im 
DE 
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矩阵 M 的 右 和 左 特征 子 空间 , Sah X YAA 
量 张 成 , 归 一 化 后 满足 YX=1; 与 特征 子 空间 对 应 ， 
满足 式 () 的 工 和 矩阵 称 之 为 矩阵 M 的 “广义 特征 矩 
阵 ”， 其中， 广义 特征 矩阵 的 特征 值 也 是 矩阵 M 的 
特征 值 。 


T 


MX = XL, Y"M = LY” 

由 上 式 (D) 可 知 , L 可 由 下 式 求解 得 到 : 
L=Y"MX 
SCARE IE FE BE Ae RE IE T E Vi] a) ll RP IE EL A 
向 量 的 推广 , 满足 : 若 (v, 力 是 广义 特征 矩阵 工 
的 特征 对 ( 式 中 ”为 工 的 右 特 征 向 量 )， 则 (CXy， 力 为 
矩阵 M 的 特征 对 。 其 中 ， 特 征 值 是 维 数 为 1 的 特 
殊 广 义 特征 和 矩阵， 特征 向 量 张 成 的 空间 是 特征 子 
空间 的 特例 ， 对 应 L 为 对 角 和 矩阵 。 为 方便 理解 及 

后 文 推导 ， 下 文 举例 说 明 。 

假设 矩阵 MM 有 一 组 特征 子 空间 R(X,) 、R(Y) 
(二 1,.…,n)， 对 应 有 nn 个 2 维 广义 特征 矩阵 Li， 表达 
式 如 式 (3) 所 示 ， s X=(Xa, X2) Y=(Yn, Yn)» 
Xas Xos Yas YER?” ., 根据 特征 子 空间 的 定义 
Hy A, AF RAM... 如) 和 (了 Pi... 也) 可 将 M 转换 
为 块 对 角 和 矩阵 ， 且 对 角 块 为 广义 特征 矩阵 L;， 表 示 
为 : 


(1) 


(2) 


特征 


(Y,...Y,)" M(X....X,,)=diag(L, ) 

RP: ¥"MX, =L 
由 于 (Xi ... XY) 和 (了 ... Yn) AEE Fe E AE 
且 列 向 量 分 别 线 性 无 关 , 故 式 (3) 的 矩阵 变换 不 改变 
矩阵 M WHEE, BIERE diag( L) 35 M 的 特征 值 
完全 相同 ， 其 中 L EIAN M 的 部 分 特征 值 ， 
Li 可 理解 为 M 经 过 特征 分 解 得 到 的 款 个 解 耦 子 矩 
阵 。 特 别 地 ， 若 M 有 nn 个 2 重 特征 值 4i=1,...,7) 
且 存 在 两 个 线性 无 关 的 特征 向 量 ， 那 么 其 张 成 的 空 
间 为 特征 子 空 间 ， 此 时 广义 特征 矩阵 满足 Liho 
1.2 ”特征 子 空间 摄 动 

考虑 如 下 受 扰 和 矩阵: 

M=M+E 

Hp, ERE E ERE M 所 受 扰动 。 
若 扰动 矩阵 E EED, MARIETTE M 
的 特征 子 空间 存在 近似 关系 ， 如 下 引 理 。 

引 理 120 (特征 子 空间 摄 动 ): WEE M 的 右 
和 左 特征 子 空间 左右 分 别 为 (对) 和 R(Y)， 那 么 
m RETRE R(X) A RY) AR SUREE 
了 的 近似 满足 : 

R(X) R(X +0(||E])). 3(Z)=sR(Z+ollEl) © 

L=Y" MX =Y" MIX +ol(|E| ) (6) 


(3) 


(4) 


# 


期 刘 晨 曦 等 ， 基 于 特征 子 


中 , ol.) 表 示 高 阶 无 穷 小 , EREE 
由 引 理 1 可 知 , M 的 特征 子 空间 可 以 近似 站 的 


Nit 


特征 子 空 间 。 扰动 越 小 , 它们 的 特征 子 空间 越 接近 。 
当下 的 特征 子 空间 计算 非常 复杂 时 ， 可 构造 合适 的 


M 简化 及 的 广义 特征 和 矩阵 计算 。 
维和 矩阵 的 例子 说 明 上 述 结论 。 
1.3 ”特征 子 系统 

考虑 2n 维 的 动态 系统 A4， 记 闭环 传递 函数 矩 
阵 和 系统 矩阵 分 别 为 Ti(s) 和 用 =T(s)， 那 么 特征 方 
程 det(M )=0 决定 了 4 的 稳定 性 。 根 据 定义 1 并 
类 比 式 (3) 可 知 , 站 有 一 组 特征 子 空间 和 多 个 低 维 广 
AURERE G, ， 且 互 的 特征 函数 也 是 让 的 特征 


附录 A 给 出 了 两 
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首先 ， 考 虑 所 有 线路 阻抗 比 相同 ， 从 装备 端口 
往 交 流 电 网 看 ， 在 全 局 同步 坐标 系 下 的 电网 端口 电 
压 微 增 量 AUiy(s) 和 电流 微 增 量 ALwy(s) 的 灵敏 度 方程 
可 表示 为 L151: 


-Ar =(B®E(s))AU, (8) 


Bimmer) 
Ve+l ke +t) 1 | 
(a+r) til -1 s/(@ +7) 
其 中 ,，B KEEMA T BS CS DB R BE 


B 


(ntm+r)l Ea 


数 能 够 决定 系统 在 某 一 关注 模 态 下 的 稳定 性 。 
利用 低 维 的 天 分 析 丰 矩阵 的 性 质 ， 引 入 了 特 生 
系统 的 概念 ， 有 具体 为 ; 
定义 2 (系统 卸 阵 和 特征 子 系统 ): 设 动态 系统 
4， 其 系统 闭环 传递 函数 逆 和 矩阵 称 之 为 该 系统 的 系 
SEE, Bl w=T(s); 记 让 的 第 大 个 特征 子 空间 和 广 
义 特征 矩阵 分 别 为 只 (总 ) RK) ALL, GER 
(M, RREH L 的 低 维 子 系统 称 之 为 4 的 特 
征 子 系 统 。 
L, =Y" MŠ „k ef{l...,n) (7) 


2 系统 建 模 及 问题 描述 
设备 类 型 多 样 运行 方式 多 变 


函 
数 。 因 此 , 广义 特征 抢 阵 五 (hE {L… 妆 ) 的 特征 函 
为 
F 子 


px 电厂 


FRE; r=R/X 为 线路 阻抗 比 ; 表示 变量 的 增 量 ; Q 
表示 Kronecker 积 ; wo 为 交流 电网 的 同步 转速 。 当 
网 络 侧 包 含 电 力 电子 无 功 补偿 、 构 网 型 装备 或 负 蓓 
等 元 件 时 ， 修 改 交 流 网 络 端口 电压 对 电流 的 灵敏 度 
矩阵 即 可 考虑 这 些 元 件 带 来 的 影响 ， 后 续 强 度量 化 
方法 仍然 适用 ， 限 于 篇 幅 不 再 详细 展开 。 

其 次 ， 装备 端口 的 电压 微 增 量 AUwy(s) 和 电流 微 
增 量 ALw(s) 的 灵敏 度 方程 可 表示 为 91; 

Yam (5) 


Yine(s) 
其 中 ，Yisr1(s) ~ Yerom (ORIA ntm 台 装 备 的 
纳 传 递 函数 矩阵，0>ER2 > ， 元 素 均 为 0。 
最 后 , 结合 式 (8) 和 (9)， 利用 如 下 行列 式 可 表示 

图 1 所 示 系 统 的 闭环 特征 方程 : 
det (T,,, (5))=det (Ying (s)+Y,,(s))=0 (10) 
其 中 ，77e(s) 为 原 系 统 的 闭环 传递 函数 逆 和 矩阵 ， 即 


aq 


tte 
` TIOR 


图 1 多 类 型 新 能 源 馈 入 电力 系统 
Fig.1 The new power System multiple types of renewa- 
ble energy devices integrated 


考虑 图 1 所 示 新 型 电力 系统 ， 包 含 ntm AR 
备 及 接 入 的 母线 ,7 个 中 间 无 源 母 线 和 1 个 无 穷 大 
母线 ( 若 有 构 网 装备 , 无 穷 大 母线 可 以 没有 ), 其 中 ， 
n 台 装 备 是 新 能 源 或 储 能 的 跟 网 型 变 流 器 ; m 台 是 
构 网 型 装备 ， 如 构 网 型 风机 、 储 能 或 调 相 机 等 。 
2.1 系统 动态 建 模 

在 全 局 xy 坐标 系 下 ， 建 立 图 1 所 示 系 统 的 线 
性 化 模型 ， 得 到 系统 的 闭环 特征 方程 。 


原 系统 矩阵 ，det(.) 表 示 求 行列 式 。 
2.2 ”基于 广义 短路 比 的 系统 强度 量化 思路 回顾 

式 (10) 所 示 闭 环 特征 方程 决定 了 系统 稳定 性 ， 
然而 ， 由 于 Tseds) 维 数 较 高 ， 故 寻找 近似 的 低 维系 
统 以 简化 分 析 。 

WMTW AR, HAA RSMAS SER 
单 馈 入 系统 。 基 于 失 稳 风险 最 高 的 等 效 单 馈 入 系统 ， 
类 比 传统 短路 比 定义 ， 可 给 出 原 系统 的 广义 短路 比 
gSCR 与 装备 / 场 站 临界 短路 比 SCRo, 并 形成 基于 广 
义 短路 比 的 系统 强度 分 析 方 法 ， 而 且 上 共有 电网 动态 
与 装备 动态 分 离 的 特点 。 计 算 公式 如 下 式 所 示 品 

SSCR> CgSCR, = SCR, aD 
g a 8SCR-SCR, 
SCR, 


或 B%* 100% > By % (12) 
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FEA, CgSCRo 为 广义 短路 比 临界 值 ; fo% 是 为 了 保 
证 系统 具有 足够 的 裕 度 而 设置 的 阔 值 ，gSCR 的 表 
达 式 为 : 


aan zl 1 E 
BSCR È Arin Sint S pbs) ] =o" [人 


式 中 : Sp0+m =diag(SB1. Seorm)，diag() 表 示 对 角 拢 
BE, Sei., Sgtmw 分 别 为 ntm 台 装 备 的 额定 容量 ; 
Brea 为 保留 装备 母线 的 网 络 导 纳 矩 阵 或 电压 电流 灵 
DUZER: Zed Bredt: F(+) All anO 分 别 表示 求解 
矩阵 的 最 大 奇异 值 和 最 小 特征 值 。 

广义 短路 比 是 灵敏 度 和 矩阵 的 最 大 奇异 值 ， 反 映 
了 多 端口 交流 网 络 中 电压 扰动 对 电流 扰动 的 最 大 
REM, 
对 于 异 构 系统 ， 利 用 和 矩阵 摄 动 理论 对 gSCR 和 
SCRo 进行 近似 即 可 。 具体 地 , 根据 摄 动 来 源 以 及 装 
备 和 电网 动态 特性 分 为 两 类 : 1) 跟 网 装备 动态 : 跟 
网 型 变 流 器 外 特性 存在 差异 ， 认 为 对 同 构 系 统 装备 
动态 存在 摄 动 ， 此 时 计算 gSCR 的 边界 条 件 不 变 ， 
而 SCRo 取决 于 装备 动态 特性 的 加 权 平 均 ， 有 具体 见 
文献 [7][17][19][20];， 2) 等 效 交 流 网 络 动态 : 当 考 
虑 构 网 型 变 流 器 或 调 相 机 等 电压 源 特 性 占 优 的 接 
入 装备 ， 可 认为 是 对 网 络 修正 ， 此 时 SCR Wits 
方法 不 变 ， 而 修正 gSCR 的 计算 边界 条 件 即 可 计 及 
这 些 构 网 装备 的 影响 ， 有 具体 见 文 献 [9][10][11]。 
2.3 拟 解决 的 问题 

综 上 所 述 ，gSCR 指标 及 其 临界 值 的 定义 和 计 
算 ， 其 关键 步骤 是 寻找 等 效 的 低 维 单 馈 入 系统 ， 即 
特征 子 系统 。 然 而 ， 现 有 研究 仅 适 用 于 单一 典型 场 
景 下 的 问题 ， 如 仅 考 虑 外 特性 异 构 的 跟 网 型 变 流 器 
I9、 构 网 型 装备 00001] 或 储 能 系统 09 等 装备 接 入 的 场 
景 ， 其 方法 及 理论 缺乏 普 适 性 ， 导 致 对 广义 短路 比 
的 理解 相对 复杂 。 为 此 ， 本 文 利 用 更 具 物 理 意义 的 
特征 子 系统 概念 ， 统 一 广义 短路 比 的 导出 原理 及 计 
算 方法 ， 主 要 讨论 两 个 问题 : 

问题 1， 近 似 原 大 规模 系统 模 态 的 低 维 解 耦 子 
系统 〈 即 特征 子 系统 ) 的 理论 基础 。 

问题 2， 基 于 特征 子 系统 ， 广 义 短路 比 及 其 临 
界 值 的 统一 导出 原理 及 计算 方法 。 
3 ”特征 子 系统 及 其 近似 方法 

结合 定义 2 可 知 , MA EMA AAR KAI 
RE, 存在 低 维特 征 子 系统 与 广义 特征 矩阵 。 其 中 ， 
特征 子 系统 可 近似 原 动 态 系统 在 关注 模 态 下 的 稳 
定性 ， 且 维 数 更 低 ， 分 析 更 简单 。 

上 述 分 析 思 路 的 关键 是 获得 特征 子 系统 , 即 求 
解 广义 特征 和 矩阵。 为 此 , 本 节 首 先 利用 和 矩阵 的 舒 尔 
补 变换 简化 原 系统 矩阵 ; 其 次 , 基于 引 理 1 获得 系 


(13) 


in 
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统 矩 阵 的 近似 特征 子 空间 , 并 得 到 对 应 的 广义 特征 
矩阵 及 低 维 特征 子 系统 。 
3.1 系统 矩阵 及 其 简化 

考虑 图 1 所 示 电 力 系统 , 下 文 从 小 扰动 同步 视 
角 给 出 其 特征 子 系统 的 近似 求解 方法 ， 其 余 横 态 下 
的 特征 子 系统 求解 方法 类 似 。 

将 系统 看 成 由 两 部 分 组 成 , 一 部 分 为 等 效 网 络 ， 
包含 交流 网 络 与 m 台 构 网 装备 的 等 效 导 纳 ( 在 分 析 
频段 内 的 外 特性 由 戴 维 南 等 效 电路 主导 ， 其 等 效 导 
纳 为 Biaw (=1,.…,m))， 另 一 部 分 包含 n 台 跟 网 装 
备 动态 模型 。 式 (10) 所 示 的 原 系统 矩阵 可 变换 为 : 

T (5)=Yine (8) +¥ur (s) (14) 
AH: YO (s)=(B+ Bx) ®6(s) 


ui 


Vise (s) 
O (p) ve 
Yi (s) 7 Y rn (s) 
Oimn) 
0 
aap 1) T, (s) T, (s) 
Ber = diag (Ba) ,Da (s)= T (s) T (s) 
0, 3 4 

其 中 ， Bsr 为 m AWERI E R BLP FB E 


的 修正 ， Tilo TAs TOM Ta(s) 为 具有 合适 维 

FIERE: DJ 2 维 单位 矩阵 ， On 0; 分别 为 n、 

r AFERE. 
假设 未 接 装 备 时 电网 是 稳定 的 ， 则 根据 Schur 

we FEA, det(Tne(s))=0 的 零点 由 det(Ti(s)- To(s) Ty 

\(s)T3(s)) =0 决定 ， 故 det(T?,.(s))=0 决定 了 图 1 

所 示 系 统 的 稳定 性 ， 其 中 (推导 见 附录 B): 
det(T, (s)—T, (s)T,"(s)T,(s)) = det (TO (s)) 
= det (YQ (s)+Y® (s)) =0 


~ (s ) l 


(5) 


式 中 ， YA(s)= : | 
Yin (8) 
YO (8) = Bia 86 (s) 
ik B, 
B+Bip= 
B, B, 
其 中 ， BIER””, BERT), BERO, B, 
ER ONM) 。 
进一步 考虑 运行 点 的 影响 ,将 剩余 n 台 跟 网 型 
装备 实际 工 况 与 其 动态 解 看 ?9。 具 体 地 ， 用 包含 装 
备 实际 工 况 信息 的 矩阵 (U2P)@Jb AH TOnei(s) 
TQ%4als) 可 转化 为 如 下 所 示 的 Tels) $ M= 
Tres), M 的 特征 函数 决定 了 图 1 系统 的 稳定 性 。 


|e =(B,-B,B;'B,) 


第 URES: ETHET 
aet(((U7/P 3 n (5))=dee(722(s)) 
=det( YG (s)+Y0 (s))=0 
式 中 : ¥O(s)= es @¢(s) 
¥{2(s) =((v2/P)@1,)¥22(s)= 人 全 Ypa oj 
其 中 ，U3/P= diag(U72/Pi)，Ui、Pi; 分 别 为 装备 i 的 端 


(16) 


有 功 功率 为 负 。 
JE 1: 当 装 备 在 额定 运行 点 , 其 端口 电压 为 1.0 
pu 输出 有 功 功 率 等 于 其 额定 容量 Si (i=1,.….,n)， 


此 时 [2/P= diag(U2/P)-diag(Sw)=Sw1， 故 额定 运 
行 点 的 分 析 方 法 是 实际 工 况 下 的 特例 。 
3.2 ”特征 子 系统 近似 方法 

on 台 等 效 装备 动态 模型 同 构 时 ， 等 效 导 纳 矩 
ME(O?/P)B rea 9h KF ARB eee 
HAM 的 一 组 特征 子 空间 。 其 中 ， 当 系统 中 有 
(ga) 台 装备 有 功 功率 反 向 (为 负 值 ) 时 ,(LP/P Be 
OLA q 个 负 特 征 值 09， 将 特征 值 排序 为 14.j<.…<4- 
OUAS.. Ang 21 为 最 小 正 特征 值 ; R(X) , 
R(V,) G= -q,...-1,1,.. 0- (U7IP)B rea Qo HHE 
值 力 对 应 的 特征 子 空间 。 进 一 步 ， 根 据 定 义 1 和 2 
并 类 比 式 G3)， 利 用 R(X,) M R(Y,) (kE 人 9 
1,1,.….,n-q}) 可 求解 原 系统 的 多 个 特征 子 系统 ， 分 别 
近似 原 系 统 不 同 模 态 下 的 稳定 性 ， 总 结 为 如 下 引 理 


口 电压 、 输 出 有 功 功 率 , 当 考 虑 装备 的 功率 反 向 时 ， 
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G,(s)=Y, ‘| diag ((U? /B) Vinay (s )) | 
=) pu (U? /P)Y pri (s), Pri T Xi Yri 
其 中 ， Thom 为 对 应 Ar SRF RY AEE; Gr 
(5) 为 等 效 加 权 装 备 模 型 ;pj 为 第 i 台 装 备 参与 因子 ， 
Xkix Yri DBA Xk~x Vk 的 第 i | 个 元 素 。 

其 次 ， 观 察 构造 的 同 构 系 统 Mi 中 ， 装 备 动态 
是 原 系统 中 所 有 接 入 跟 网 型 装备 动态 的 加 权 。 再 根 
据 引 理 1 和 2 TY, R(X,) 和 网 (又 ) 仍然 是 M: 的 
特征 子 空间 ， 利 用 该 组 特征 子 空间 可 求解 政 的 广 
SURES ML, ， 表 示 为 : 

L, ~Y MX, =G,(s)+4,-¢(s) (19) 
A,-$(s)=¥," ((U?/P) Bi. Ge 总 
ye (U°/P)B ax, (5) 

最 后 , 将 五 看 作 一 低 维 系统 的 系统 矩阵 , 那么 
该 低 维 系统 是 原 系统 的 一 个 特征 子 系统 , 可 近似 原 
系统 在 对 应 模 态 下 的 稳定 性 (由 于 上 上 ,保留 了 原 系 
统 和 矩阵 的 零 极 点 ， 故 还 能 反映 原 系统 的 鲁 棒 性 )。 
该 结论 总 结 为 定理 1， 具体 证 明 见 附录 。 

定理 1: 将 原 系统 矩阵 让 看 作 同 构 系 统 矩 阵 
Mi 的 摄 动 ，Mi 的 特征 子 空间 与 让 的 特征 子 空间 近 
似 相 等 ， 且 六 的 广义 特征 矩阵 忆 可 由 式 (19) 计 算 。 
上 述 得 到 的 特征 子 系统 维 数 较 低 , 并 能 够 近似 


A: 


2， 其 具体 证 明 见 附录 C。 
引 理 2: 广义 特征 矩阵 五 满足 : 
L =Y" MX, =Y pra (5 s)+A,-¢(s) (17) 
RP: M =I, 8Y pr (8) +(U?/P)Bi,, @6(s) 
其 中 , Yigres(s) 为 相似 的 单 台 等 效 装备 动态 模型 ;Xi= 
“Oh, Y= OR, Xi. Yh WW ea HA 
(W7/P)B redQb RF AWA. ARES, xi ye E 
R”, x ye 分 别 为 (UWP)B rea KF AWA EFF 
征 向 量 。 
当 装 备 动态 模型 存在 差异 时 ,构造 等 效 同 构 系 
统 , 其 系统 窍 阵 的 特征 子 空间 可 近似 原 系统 矩阵 的 
特征 子 空间 ， 再 由 引 理 1 可 得 广义 特征 矩阵 五 的 
HEAL. Alpha, HERH R(X, ) F RY, ) 构造 等 效 
同 构 系统 矩阵 Mi 为 : 
M, =D (s)=Y r (5)+ Yer hom (s) (18) 
RP: Yoo jum (8)=(U?/P) Bi, 86 (5) 


Vier (s), OG, (s) 


反映 原 系 统 在 对 应 模 态 下 的 特性 。 因 此 ， 基 于 特征 
子 系统 简化 了 复杂 系统 的 分 析 。 

注 2: 当 跟 网 装备 动态 模型 同 构 ， 等 效 同 构 系 
统 矩 阵 与 原 系统 矩阵 的 特征 子 空间 完全 相同 ， 特 征 
子 系统 稳定 性 与 原 系统 相同 ;， 反 之， 等 效 同 构 系统 
矩阵 与 原 系统 和 矩阵 间 存 在 差异 , 广义 特征 矩阵 的 近 
似 与 真实 广义 特征 矩阵 存在 偏差 此 时 特征 子 系统 
近似 反映 了 系统 稳定 性 和 和 鲁 棒 性 。 
3.3 ”特征 子 系统 的 物理 意义 
据 上 节 可 知 ,利用 特征 子 空间 摄 动 可 求解 原 
系统 矩阵 的 广义 特征 矩阵， 对 应 的 特征 子 系统 可 
用 于 分 析 原 系统 稳定 性 ; 物理 上 看 , 特征 子 系统 是 
一 个 包含 装备 与 网 络 动 态 的 低 维 单 馈 入 系统 ， 本 
质 是 原 系统 在 模 态 坐标 系 下 的 一 个 解 耦 子 系统 。 特 
征 子 系统 的 物理 意义 及 求解 过 程 可 用 图 2 和 图 3 
总 结 
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由 式 (19) 与 图 2 可 知 , 特征 子 系统 为 单机 无 穷 
大 系统 ， 由 跟 网 型 等 效 加 权 装 备 与 等 效 线路 构成 。 
其 中 ， 接 入 装备 的 特性 为 原 系统 中 跟 网 型 装备 的 加 
权 平 均 ， 反 映 系统 中 装备 动态 的 整体 特性 ， 等 效 线 
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路 动态 计 及 了 构 网 装备 等 效 导 纳 等 因素 对 网 络 动 
态 的 影响 。 进 一 步 ， 从 传统 短路 比 看 ， 特 征 子 系统 
的 短路 比 取 值 为 对 应 等 效 线路 导 纳 ， 其 大 小 为 特征 
(E Ar 的 绝对 值 ， 可 用 于 表征 特征 子 系统 的 电网 强 
度 。 


实际 复杂 系统 


跟 网 型 装备 1 


特征 子 空间 [| Ó 稳定 性 
近似 。 “近似 相同 


第 1 个 关注 模 态 下 的 特征 子 系统 


第 1 个 跟 网 型 
等 效 加 权 装 备 。 等 效 线路 


ot 


G(s) 


SCR= |4 | 


第 个 关注 模 态 下 的 特征 子 系统 


第 n 个 跟 网 型 等 效 线路 


等 效 加 权 装 备 
-E © 
SCR = |. | ! 压 源 


Gs) 


子 系统 物理 意义 示意 图 


图 2 特征 


Fig.2 The schematic diagram of physical significance of 
eigen-subsystem. 


基于 特征 子 空间 摄 动 的 特征 子 系统 求解 


Treks) 
见 式 (8)-(10) 


获得 原 系统 物理 模型 


和 矩阵 拆 分 、Schur 补 
得 到 简化 系统 矩阵 用 
见 式 (14)-(16) 


根据 稳定 模 态 ， 简 化 跟 
网 或 构 网 装备 动态 


近似 模 态 解 耦 


RETAS: WERK 
特征 子 系统 矩阵 L ( 见 式 (19)) 
注 : L, 为 广义 特征 和 矩阵 近似 


基于 特征 
求解 特征 子 


F 子 空间 摄 动 
空间 近似 


图 3 特征 子 系统 求解 方法 示意 
Fig.3 The schematic diagram of the calculation method of 
eigen-subsystem. 


还 值得 注意 的 是 ， 由 于 不 同 模 态 的 稳定 机 理 不 
同 ， 其 对 应 的 特征 子 空间 及 特征 子 系统 应 不 同 ， 对 
应 下 取 值 不 同 。 因 此 ， 实 际 问题 中 需 根据 具体 模 态 
及 稳定 机 理 ， 选 择 相应 的 特征 子 系统 ， 实 现 复杂 系 
统 的 简化 分 析 。 此 外 ， 由 于 本 文 聚焦 跟 网 型 装备 的 
主导 模 态 并 将 构 网 装备 等 效 为 对 地 支 路 ， 对 应 特征 
子 系统 仅 为 跟 网 型 装备 的 单机 无 穷 大 系统 ， 然 而 ， 
在 其 余 模 态 下 ， 对 应 特征 子 系统 形式 更 为 复杂 ， 包 
括 但 不 限于 单机 无 穷 大 系统 。 


4 ”基于 特征 子 系统 的 广义 短路 比 导 出 原理 
及 强度 量化 方法 

本 节 提 出 了 基于 特征 子 系统 的 广义 短路 比 统 
一 导出 原理 ， 并 给 出 了 相应 流程 和 几 种 典型 场景 下 
的 计算 方法 。 
4.1 ”特征 子 系统 中 的 电网 强度 和 广义 短路 比 指标 

根据 3.3 节 可 知 ， 特 征 子 系统 可 近似 表征 原 系 
统 的 稳定 性 , 并 由 等 效 加 权 装 备 动态 与 等 效 线路 动 
态 构成 ， 电 网 短路 比 大 小 为 特征 值 丰 的 绝对 值 。 因 


此 ， 原 系统 强度 取决 于 特征 子 系统 动态 ， 当 等 效 加 
权 装 备 动态 一 定时 ， 原 系统 强度 取决 于 特征 子 系统 


的 短路 比 ， 即 等 效 线路 导 纳 。 

当 电网 电压 支撑 强度 低 时 ， 静 态 电压 稳定 及 由 
锁 相 环 带 来 的 小 扰动 失 稳 风险 较 高 ， 故 这 两 类 稳定 
问题 所 对 应 的 特征 子 系统 ( 即 所 1) 备 受 关注 。 具 
体 地 ， 在 该 模 态 下 ， 短 路 比 越 小 ， 系 统 失 稳 风险 越 
高 中; 那么 对 应 短路 比 为 如 (对 应 厂 1) 的 特征 子 系 统 
最 易 发 生 失 稳 , 该 短路 比 被 定义 为 “广义 短路 比 ”， 
如 式 (20) 所 示 。 因 此 ， 广 义 短 路 比 仅 是 基于 特征 二 
系统 的 一 个 具体 应 用 ， 其 余 重 要 的 特征 子 系统 ， 也 
能 够 反映 原 系统 在 其 余 模 态 下 的 特性 ， 而 且 都 能 提 
出 对 应 的 指标 。 然 而 ， 般 来 说 ， 基 于 短路 比 的 分 
析 方 法 适用 于 静 稳 及 小 扰动 同步 等 模 态 对 应 的 特 
征 子 系统 ， 其 特点 是 这 类 问题 与 电网 电气 距离 强 相 
关 。 至 于 稳定 机 理 为 LC 谐振 的 高 频 振荡 对 应 的 特 
征 子 系统 ， 基 于 短路 比 的 分 析 方 法 可 能 不 再 适用 ， 
具体 见 文献 [7] 的 讨论 。 


gSCR=y! (U’/P) B aX in (U °/P) B'a) 
= Aig (UP) Be (U"/P | 
=o" G |p)? Zu (U? ie)" | 


其 中 ， AK, (-) ea READ TERRIA, Ze 


> 


(20) 


B red! o 
此 外 , 由 式 (20) 可 知 , J AEK EE BE BY LURTE 


当 闭 环 系统 的 特征 子 空间 与 开 环 系统 〈 即 仅 考虑 电 
Px) SRE RE) 的 特征 子 空间 相同 时 , 如 同 构 系统 中 ， 
两 者 结论 一 致 ， 反之 ， 基 于 特征 值 的 定义 仍然 是 一 
种 误差 较 小 且 便 于 实际 操作 的 计算 方法 。 
再 者 , 结合 式 (19) 和 (20) 还 可 知 ， 当 等 效 加 权 装 
备 / 场 站 动态 一 定时 , 根据 广义 短路 比 与 系统 强度 的 
正 相 关 性 , 能 够 确定 等 效 加 权 装 备 / 场 站 在 达到 临界 
稳定 状态 或 给 定性 能 要 求 所 需 的 最 小 短路 比 (简称 
“等 效 加 权 装 备 / 场 站 临界 短路 比 SCRo”)。SCRo 可 
反映 装备 / 场 站 动态 对 电网 强度 的 耐 受 能 力 , 恰恰 就 
是 广义 短路 比 临 界 值 CgSCRo ， 表 示 为 : 

SCR, = arg {det (G, (s,)+gSCR-G(s,))=0} (21) 


其 中 ，arg{} 表 示 求 方程 的 根 ; su 为 等 效 加 权 装 备 在 
给 定性 能 要 求 下 解 得 系统 主导 特征 根 。 

综 上 ， 复 杂 系 统 的 广义 短路 比 统一 导出 过 程 为 : 
1) 利用 等 效 同 构 系 统 特征 子 空间 ， 求 解 原 系 统 矩 
BE TEE ALLA SCRE IE FE ME, 获得 能 够 近似 原 系 统 动态 
特性 的 特征 子 系统 ; 2) 基于 特定 的 特征 子 系统 , 计 
算 其 短路 比 指标 ( 即 广 义 短 路 比 ) 及 装备 / 场 站 临界 
值 ( 等 于 广义 短路 比 临 界 值 ), 进而 形成 融合 电网 ) 
义 短 路 比 与 装备 / 场 站 临界 短路 比 的 量化 方法 。 其 中 ， 
广义 短路 比 仅 与 等 效 网 络 动态 和 接 入 装备 的 容量 
有 关 、 临 界 值 仅 与 等 效 加 权 装 备 / 场 站 的 动态 有 关 ， 
故 广 义 短路 比方 法 保留 了 传统 短路 比方 法 所 具有 
的 网 络 动态 与 装备 动态 分 离 量化 的 优势 。 
4.2 ”典型 场景 下 广义 短路 比 计算 方法 

结合 3.1 节 中 系统 矩阵 的 简化 方法 、 公 式 ( 见 式 
(15)-(16)) 以 及 4.1 节 中 广义 短路 比 的 计算 公式 〈 见 
(20)), 本 节 给 出 三 种 典型 场景 中 系统 和 矩阵 的 具体 形 
式 ， 由 此 实现 广义 短路 比 的 计算 。 
4.2.1 考虑 装备 在 非 额定 运行 点 的 广义 运行 短路 比 

考虑 图 1 所 示 系 统 仅 包含 2 台 采 用 不 同 控制 策 
略 、 参 数 的 新 能 源 变 流 器 ，m=0。 其 中 ， 新 能 源 变 
流 器 在 不 同 的 非 额定 运行 点 ， 输 出 有 功 功 率 及 端口 
电压 分 别 为 0.2~1p.u 和 0.9~1.1p.u. 不 等 231。 根 据 式 
(15)-(16), (20), 即 可 求 得 该 场景 下 简化 后 的 系统 算 
阵 及 广义 短路 比 。 具 体 过 程 为 : 

首先 ， 考虑 无 构 网 装备 等 效 导 纳 对 网 络 动 态 的 
修正 , 即 Bzsr=0w+i, 将 Bisr=0n4r 代 入 式 (15) 并 令 m=0， 
可 求 得 仅 包 含 n 台 装备 接 入 母线 的 网 络 等 效 导 纳 算 
阵 Breas HK, 考虑 将 n 台 装 备 的 实际 工 况 与 其 动 
SAR YO (s) 与 U*”/P 如 式 (22) 所 示 ; 最 终 ， 将 


阵 特 征 值 定义 ， 也 能 通过 闭环 系统 的 特征 向 量 得 到 。 
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Breas YO (s) 与 UWP 分 别 代入 式 (16)、(20)， 可 求 
得 该 场景 下 的 系统 矩阵 与 广义 短路 比 。 

Yige (8) = diag (Ysa (8)),U?/P = diag (U? /P) (22) 
其 中 ，Yysci(s) 为 第 i 台新 能 源 变 流 器 导 纳 传递 函数 
ERE, i=1,.….,n; Un Pi 分 别 为 装备 i 的 端口 电压 、 
输出 有 功 功率 。 
4.2.2 考虑 装备 有 功 功率 反 向 的 广义 短路 比 

根据 装备 的 物理 特性 与 系统 需求 ， 部 分 装备 可 
反 向 从 电网 吸收 功率 ， 如 储 能 系统 与 柔性 直流 输电 
等 。 不 失 一 般 性 地 ， 以 常规 储 能 系统 为 例 ， 考 虑 图 
1 所 示 系 统 包含 n 台 采 用 不 同 控制 策略 、 参 数 的 新 
能 源 变 流 器 ，m 台 常 规 储 能 系统 。 其 中 ， 装 备 均 运 
行 在 额定 工 况 ; 新 能 源 变 流 器 主要 向 电网 发 出 有 功 
功率 ， 储 能 变 流 器 主要 向 电网 发 出 (放电 ) 或 吸收 有 
TH (充电 )。 根据 式 (15)-(16)、(20), 即 可 求 得 该 场景 
下 简化 后 的 系统 矩阵 及 广义 短路 比 。 有 具体 过 程 为 : 
首先 ， 考 虑 无 构 网 装备 等 效 导 纳 对 网 络 动态 的 
修正 ， 即 Bisr=0nimir:， 将 Bir =nimi 代入 式 (15)， 
可 求 得 仅 包含 ntm 台 装 备 接 入 母线 的 网 络 等 效 导 
内 矩阵 Breas 其次， 考虑 将 ntm 台 装 备 的 额定 容 
量 与 其 动态 解 厢 ， 其 中 ， 当 储 能 变 流 器 充电， 其 额 
定 容量 需 加 上 负 号 , PIM, 考虑 m 台 储 能 均 工作 在 
RERA YO (s) 5U /P WRITER: RA, 将 
B'as Yigg (8) 与 UP 分 别 代 入 式 (16)、(20)， 可 求 
得 该 场景 下 的 系统 矩阵 与 广义 短路 比 。 

YƏ (s) z ~ (s)) 


五 


ALS 


lin 


diag (Yocy w) (23) 


diag (sxi) 
—diag (Sci) 


其 中 ，Ypcg(s) 为 第 j 台 储 能 变 流 器 的 导 纳 传递 函数 
HERE, j=1,... m; Swix Smt 分 别 为 第 i 台新 能 源 
变 流 器 及 第 7 台 储 能 变 流 器 的 额定 容量 。 
4.2.3 考虑 构 网 型 装备 影响 的 广义 短路 比 

考虑 图 1 所 示 系 统 包 含 n 台 采 用 不 同 控制 策 
略 、 参 数 的 新 能 源 变 流 器 ，zz 台 构 网 型 储 能 系统 或 
调 相 机 。 其 中 ， 装 备 均 运 行 在 额定 工 况 ; 调 相机 和 
外 特性 近似 为 电压 源 且 惯量 系数 较 大 的 构 网 型 变 
流 器 能 够 为 系统 提供 电压 支撑 能 力 ， 在 锁 相 环 主导 
的 小 扰动 同步 稳定 问题 中 能 够 等 效 为 含 内 电抗 的 
电压 源 。 根据 式 (15)-(16)、(20)， 即 可 求 得 该 场景 下 
简化 后 的 系统 矩阵 及 广义 短路 比 。 上 有 具体 过 程 为 : 
首先 , 考虑 m 台 构 网 装备 等 效 导 纳 对 网 络 动态 
的 修正 ， 即 Bigr= diag(0 diag(Ber) 0 ; )， 将 其 代入 
式 (13)， 可 求 得 仅 包含 n 台 装 备 接 入 母线 的 网 络 等 
效 导 纳 矩 阵 B’eq; 其 次 ,考虑 将 n 台 装 备 的 额定 容 


U’/P > Sarm) = 


Co 


中 国 电 


EIAM YO (s) 与 UV*/P 如 式 (24) 所 示 ; 最 
,将 Breas YO (8) 与 U0?/P 分 别 代 入 式 (16)、(20)， 
可 求 得 该 场景 下 的 系统 矩阵 与 广义 短路 比 。 
Was )= diag ( Ysa (s )) 
U?/P=Sy, =diag (Sy, i) 


tP, Bey HHE j 台 构 网 型 变 流 器 或 调 相 机 的 等 
导 纳 ， J=1,... 


a 
se dI 


\NS 
R ffin] 


(24) 


效 
sm, Bipr=diag(SvsmjBp 或 1/xa”), xaj” 


N 


机 


C 程 报 


为 第 7 台 调 相机 的 等 效 内 电抗 ， 即 次 暂 态 电抗 ， 
1/SvswBp 为 第 j 台 构 网 型 变 流 器 的 等 效 内 电抗 ， 
Sys 为 第 j 台 构 网 型 变 流 器 的 额定 容量 , By 为 单位 
容量 下 构 网 型 变 流 器 的 等 效 导 纳 。 

为 直观 对 比 上 述 3 种 典型 场景 下 原 系统 矩阵 的 
构成 ， 表 1 总 结 了 上 述 公 式 。 此 外 ， 表 1 还 给 出 了 
上 述 场景 对 系统 强度 的 作用 机 理 ， 详 细 讨 论 可 分 别 
参考 文献 [9][10][11][19][20]。 


表 1 典型 场景 下 原 系统 矩阵 构成 及 其 对 系统 强度 的 作用 机 理 


Tab.1 The elements required for calculation of system strength and their physical significance under typical 


scenarios 
原 系统 矩阵 M =Y (s)+¥O(s) 
Yin (s) =((U?/P) @1,)¥ine(s) Y (s)=(U"/P)B,,, ®6(s) 
型 场景 影响 
al Bret: B+Bisg 降 阶 至 YOysg(s) 端 口 的 网 络 mam 
YO er(s) UIP 等 效 导 纳 和 矩阵 
B 构 网 装备 等 效 导 纳 Bisr 
非 额 定 运行 点 diag(Yysci(s)) diag(U7/P) | 包含 所 有 装备 On» 
有 功 功率 反 向 diag(Yysci(s) Ypcs(s)) Sym)! 及 中 间 节 点 的 Onim+r 改变 gSCR BK SCRo 
构 网 型 装备 diag( Yrsci(s)) Syn" 网 络 导 纳 矩阵 diag(0, diag(Brer) 0, ) 
= 天 稳定 外 的 其 它 稳定 问题 
43 ”系统 强 度量 化 步骤 压 稳 定 外 的 其 它 稳定 问题 。 
为 方便 实际 应 用 ， 本 节 给 出 基于 广义 短路 比 的 5 算 例 分 析 及 验证 


系统 强度 量化 步骤 ， 有 具体 如 下 : 


1) 根据 装备 物理 特性 将 装备 动态 分 别 看 作 跟 
网 型 装备 动态 或 对 网 络 动态 的 修正 ; 

2) 根据 网 络 参 数 与 装备 特点 ， 基 于 式 (15) 得 到 
等 效 网 络 动态 B+Bisa; 根据 实际 工 况 、 装 备 数 和 B+ 


Bisrg， 基 于 式 (15) 获 得 等 效 导 纳 矩阵 (U2/P)B reas 

3) 根据 式 (20),， 求解 等 效 导 纳 或 灵敏 度 矩 阵 的 
特征 值 或 奇异 值 ， 即 广义 短路 比 ; 
4) 根据 式 (21) 或 (25) 求 解 SCRo。 具 体 地 ， 前 者 
可 通过 解析 计算 或 仿真 测试 等 方式 获得 门 ， 得 到 的 
SCRo 精度 较 高 ; 后 者 则 需 结 合 实际 工程 经 验 , 利用 
参与 因子 ( 见 式 (18)) 计 算 获 得 。 


SCR, =" p,SCR, , 
第 ; 台 装备 达到 临界 稳定 所 需 的 最 


(25) 


RP, SCRo ;为 
小 短路 比 。 

5) 根据 式 (12) 求 解 6%， 量 化 强度 。 

值得 注意 的 是 ， 当 系统 为 黑箱 模型 时 ， 等 效 网 
络 动 态 与 实际 工 况 矩阵 可 通过 扰动 辨识 方法 人 9 与 
PMU 量 测 等 方法 获得 数据 。 

综 上 所 述 , 通过 引入 特征 子 系统 , 统一 了 多 种 
场景 下 广义 短路 比 及 临界 值 的 导出 原理 及 计算 方 
法 ， 具 有 优势 为 : 1) 保留 物理 特性 。 特 征 子 系统 
保留 了 装备 与 网 络 的 动态 特性 ; 2) 方 法 通用 且 拓 展 
性 好 。 不 仅 适 用 于 异 构 多 馈 入 系统 的 复杂 场景 ， 而 
且 有 潜力 推广 应 用 到 除 小 扰动 同步 稳定 和 静态 电 


首先 基于 4 馈 入 系统 验证 特征 子 系统 求解 方法 
及 其 在 描述 原 系统 稳定 性 的 有 效 性 ; 其 次 ， 利 用 某 
实际 新 能 源 基 地 的 电磁 暂 态 模型 ， 验 证 广义 短路 比 
指标 在 新 型 电力 系统 中 的 有 效 性 ， 并 分 析 装 备 在 非 
额定 运行 点 、 有 功 功率 反 向 和 包含 构 网 型 装备 对 指 
标 计算 与 系统 强度 的 影响 规律 。 


5.1 ”特征 子 系统 及 其 求解 方法 的 有 效 性 

为 验证 特征 子 系统 在 描述 原 系统 稳定 性 的 有 
效 性 ， 搭 建 4 馈 入 系统 电磁 暂 态 仿真 模型 ， 拓 扑 如 
图 4 所 示 。 其 中 节点 1、4 为 常规 、 构 网 型 储 能 系 
统 接 入 节点 , 节点 2、3 为 新 能 源 变 流 器 接 入 节点 ， 


装备 控制 策略 、 参 数 及 网 络 参数 如 附 表 


E1-E3 所 示 。 


图 44 馈 入 系统 
Fig.4 The topology of a four-converter System 


人 体 地 ， 考 虑 等 比例 增 大 附 表 E3 所 示 所 有 线 
路 电感 ， 比 例 参数 k=l, 1.5. 2. 3. 4, Heb, WH 
储 能 变 流 器 从 电网 吸收 功率 ， 构 网 型 变 流 器 向 电网 
发 出 功率 , 装备 实际 工 况 如 表 2 所 示 。 以 此 为 基础 ， 
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首先 ， 分 别 基 于 特征 子 空 间 摄 动 理论 求解 特征 子 系 
统 ; 其 次 ,对比 原始 4 馈 入 系统 和 特征 子 系统 的 主 
导 特 征 根 ， 验 证 特征 子 系统 的 有 效 性 。 

首先 ， 以 算 例 1 为 例 ， 给 出 特征 子 系统 的 具体 
求解 过 程 。 有 具体 步骤 为 : 构建 系统 矩阵 ， 考 虚构 网 
型 变 流 器 对 网 络 的 修正 ， 将 系统 矩阵 降 阶 ， 并 用 


(UP) Qh 左 乘 和 矩阵 ， 最 终 系统 矩阵 变换 为 : 
det (7) (s))=det(¥2(s)+¥Q)(s))=0 26) 
-55 14 14] 
式 中 : (U?/P)B' =| -2 7 2 
4 4 16| 


3 nol Ut U; U; 
Y9 (s) =diag (Ers (s) po (s) pe oj 
APF, Frsci (8)、Yyscz (3)、Ypcs(s) 分 别 为 新 能 源 变 
流 器 与 常规 储 能 变 流 器 的 导 纳 传递 函数 矩阵 ;U2/P 
=diag(U1?/P U2?/P2 Us/P;)。 
FE RE(U?/P) Brea @ A = — ERIH 55.19, 
22=17.10, 41=-53.62, 特征 子 空间 R(X,)、R(Y) g= 
-1,1,2) 满 足下 式 , 其 中 =x 了 ,Y=yj@ h, x, 
沪 分 别 为 (U2P) Ba 关于 为 的 右 、 左 特征 向 量 。 
Y” (UI/P)B',®L)X=L BL 07) 
将 系统 矩阵 解 看 为 三 个 广义 特征 矩阵 ， 即 三 个 
特征 子 系统 的 系统 和 矩阵， 表示 为 : 
L =Y; Ta (s)X;=G,(s)+4,-¢(s) (28) 


het 


— 


式 中 : 二 
2 


U? U; U. 
G; =P; PpYs (s)+P is PY (s)+ Djs BP Ysc (s) 


3 
HF, 1,2,3, xj. vi PRA x yi 的 第 i 个 元 素 。 
表 2 算 例 的 比例 参数 及 实际 工 况 
Tab.2 The proportional parameters and actual op- 
erating conditions of examples 


> 


算 例 | k U P gSCR 
1 1 | 1.05, 0.99, 1.03 5.19 
2 1.5 1.04, 0.97, 1 3.39 
3 2 | 1.03, 0.95, 0.99 | -0.2、 1, 0.6 2.47 
4 3 | 1.01, 0.95, 0.98 1.67 
5 4 | 0.99, 0.93, 0.97 1.21 


图 5 系统 等 效 前 后 主导 特征 根 对 比 图 


Fig.5 Comparison of dominant roots before and after 


equivalence 

进一步 ， 结 合 4.1 节 分 析 可 知 ， 理 论 上 最 小 正 
特征 值 Ar 对 应 的 特征 子 系统 能 够 近似 表征 原 系统 
的 小 扰动 同步 特性 。 为 验证 该 结论 , 图 5 给 出 4 馈 
入 系统 和 11 对 应 特征 子 系统 的 主导 特征 根 。 可 以 看 
出 ， 随 着 电网 强度 的 变化 ， 原 始 系统 和 该 特征 子 系 
统 的 特征 根基 本 一 致 ， 验 证 了 该 特征 子 系统 可 近似 
原 系 统 的 小 扰动 同步 特性 。 
5.2 ”基于 广义 短路 比 的 系统 强度 量化 

搭建 菜 实际 新 能 源 基 地 的 电磁 暂 态 仿真 模型 ， 
如 图 6 所 示 。 该 系统 共 包 括 4 个 区 域 共 27 个 风力 
发 电机 组 、2 个 常规 储 能 系统 和 2 个 构 网 型 储 能 系 
统 ， 其 中 风力 发 电机 组 或 储 能 系统 的 下 标 数字 代表 
其 内 部 采用 变 流 器 编号 ; 图 中 节点 1~54 表示 装备 
接 入 节点 , 其 中 节点 43, 54 为 构 网 型 变 流 器 接 入 节 
点 ， 节 点 14、15、28、29 为 常规 储 能 变 流 器 接 入 节 
， 其 余 节 点 为 新 能 源 变 流 器 接 入 节点 ; 节点 
55~92 表示 中 间 节 点 ; 节点 93 为 场 站 外 电网 等 值 
节点 。 有 具体 地 ， 考 虑 如 下 四 个 算 例 情况 ， 算 例 的 装 
备 控制 策略 、 参 数 及 网 络 参 数 均一 致 ， 如 附 表 El- 
E3 所 示 ; 装备 实际 工 况 如 附 表 E4 所 示 。 

算 例 1: 系统 中 装备 均 运行 在 实际 工 况 ， 即 输 
出 有 功 功 率 及 端口 电压 往往 偏离 额定 值 ， 其 中 ， 常 
规 储 能 从 电网 吸收 功率 ， 其 余 装 备 向 电网 发 出 功率 ; 

算 例 2: 在 算 例 1 的 基础 上 将 装备 工 况 均 调整 
为 额定 工 况 ， 其 中 ， 常 规 储 能 仍 从 电网 吸收 功率 ， 
其 余 装 备 向 电网 发 出 功率 ; 
算 例 3: 在 算 例 1 的 基础 上 将 常规 储 能 变 流 器 
的 工作 模式 改变 为 发 出 功率 ; 

算 例 4: 在 算 例 1 的 基础 上 切除 构 网 型 变 流 器 ， 
即 系统 中 构 网 型 变 流 器 接 入 节点 悬空 。 


a 


© 
39 
x ©— 
y d Lines9.78 a 220kV 
CD 一 0.69kV 
p Pen k = 风电 机 组 
ca -on B mi 
42 
y a 构 网 型 储 能 
Lines2,78 93 


图 6 某 实 际 新 能 源 基地 拓扑 图 
Fig.6 The topology of a practical renewable energy plant 


首先 ， 按照 上 述 指标 计算 流程 分 别 解析 计算 四 
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个 算 例 中 系统 的 广义 短路 比 ， 其 次 ， 基 于 式 (21) 并 


中 


国 电 


机 


通过 搭建 等 效 加 权 单 机 仿真 或 解析 方法 得 到 四 个 


算 
gSCR 45 SCRo Ft A Ata xe 
计算 po%， 以 评估 系统 强度 。 


侈 中 系统 的 SCRo; 最 后 ， 基 于 式 (11)， 通 过 比较 
FE, 进一步 ,基于 式 (12) 


表 3 给 出 评估 结果 。 可 以 看 出 ， 算 例 1、3、4 


中 gSCR 均 大 于 SCRo， 其 中 ， 


A 


例 1 “A B%>20%, 


满足 小 扰动 稳定 裕 度 需求 ; SF 


例 2 中 gSCR=2.39 


<2.63=SCRo, 系统 失 稳 。 该 判 稳 结果 与 图 7(a)(b) 所 


示 算 例 1、2 中 装备 端口 电压 的 时 域 仿真 趋势 一 致 ， 
验证 了 基于 广义 短路 比 的 量化 方法 的 有 效 性 。 
表 3 算 例 的 指标 数值 
Tab.3 The values of indexes under four scenarios 
算 例 gSCR SCRo p% 
1 3.28 2.67 22.8% 
2 2.39 2.63 -9.1% 
3 3.03 2.68 13.1% 
4 3.02 2.69 12.3% 


gSCR=3.28>2.67=SCRo 


P i 
& 1.05 
= 
1 
0.95 — 
0 02 04 06 0.8 1 1.2 
t/s 
(a) 
12l gSCR=2.39<2.63=SCR 
A 
a 1 
D 
0.8 
0.6 moa 4 aa 
0 02 04 06 08 1 1.2 
t/s 
(b) 
L2f <a 
gSCR=2.67=SCRo 
= 1.1 
fay 
S į 
0.9 
0 0.2 04 06 08 1 1:2 
t/s 
(c) 
图 7 装备 端口 电压 的 时 域 波 形 


Fig.7 Time domain waveform of terminal voltages of 
devices in a practical renewable energy plant under differ- 
ent Operating conditions 


此 外 ,图 7(c) 给 出 算 例 1 中 系统 在 gSCR= SCRo 
时 装备 端口 电压 的 等 幅 振荡 波形 ， 验 证 了 gSCR= 
SCRo 在 描述 系统 稳定 边界 的 有 效 性 , 此 时 系统 网 络 
参数 如 附 表 E3 所 示 。 


ay 


C 程 qk 


braa 


5.3 ”广义 短路 比 的 影响 因素 
5.3.1 ， 非 额定 运行 点 的 影响 

考虑 对 比分 析 装 备 在 实际 工 况 ( 算 例 1) 与 额定 
工 况 的 场景 ( 算 例 2)。 
首先 , 在 指标 计算 过 程 中 , 算 例 1、2 的 区 别 仅 
FEF IESE UVP 的 不 同 。 具体 为 : 额定 工 况 下 UVP= 
S Ns21，SN5 为 52X52 AXA FORE, TOR AE 
定 容量 ， 当 装备 吸收 功率 ， 额 定 容 量 需 加 上 负 号 ; 
实际 工 况 下 LP/P= diag(U2/Pi)ss,，diag(U2/Pis; 为 52 
x 52 WY FA FER. 

其 次 , 如 表 3 所 示 , 算 例 1、2 中 gSCR 及 SCRo 
都 随 着 工 况 而 发 生 了 变化 ， 说 明了 实际 工 况 对 系统 
强度 的 影响 。 
5.3.2 ”有 功 功率 反 向 的 影响 

考虑 对 比分 析 常 规 储 能 变 流 器 吸收 有 功 功 率 
( 算 例 1) 与 发 出 有 功 功率 的 场景 ( 算 例 3)。 
首先 ， 根 据 常 规 储 能 变 流 器 在 充电 或 放电 工作 
模式 下 的 物理 特性 可 知 ， 在 指标 计算 过 程 中 ， 算 例 
1、3 的 主要 区 别 在 于 和 矩阵 UP 的 不 同 。 有 具体 为 ， 
常规 储 能 变 流 器 吸收 功率 时 , U2/P 中 与 其 接 入 节点 
对 应 的 对 角 元 为 负数 ， 使 得 等 效 导 纳 和 矩阵 有 对 应 的 
负 特 征 值 ， 而 发 出 功率 时 ， 相 应 对 角 元 为 正 数 ， 对 
应 特征 值 均 为 正 数 。 

其 次 , 比较 两 算 例 中 指标 及 临界 值 。 可 以 看 出 ， 
相 较 于 储 能 变 流 器 工作 于 放电 模式 ,功率 反 向 后 系 
统 的 gSCR 得 到 一 定 提升 ， SCRo 稍微 改变 , 说 明 储 
能 变 流 器 工作 于 充电 模式 会 提升 广义 短路 比 并 改 
AE SCR。。 
5.3.3 构 网 型 装备 的 影响 

考虑 构 网 型 变 流 器 接 入 ( 算 例 1) 与 不 接 入 的 场 
景 ( 算 例 4)。 
先 ， 根 据 构 网 型 装备 的 物理 特性 可 知 ， 在 指 
标 计算 过 程 中 , 算 例 1、4 的 主要 区 别 在 于 装备 对 等 
效 网 络 动态 的 修正 Bi 不 同 。 有 具体 为 ， 考 虑 构 网 型 
变 流 器 接 入 时 , Bysr 包含 构 网 型 变 流 器 的 等 效 导 纳 ， 
会 修正 网 络 动态 ， 不 接 入 时 Ber 中 元 素 均 为 0。 

其 次 ， 比 较 算 例 1、4 中 广义 短路 比 及 SCRo. 
可 以 看 出 ， 接 入 构 网 型 变 流 器 后 系统 广义 短路 比 得 
到 一 定 提升 ， SCRo 稍微 改变 , 说 明 构 网 型 变 流 器 的 
接 入 会 提升 广义 短路 比 并 改变 SCR。。 
6 结论 与 展望 

本 文 给 出 了 多 馈 入 系统 的 特征 子 系统 定义 、 物 
理 意义 及 计算 方法 ， 并 提出 复杂 场景 下 广义 短路 比 
及 其 临界 值 的 通用 导出 原理 与 计算 方法 ,包含 在 非 
额定 运行 点 、 有 功 功率 反 向 和 包含 构 网 型 装备 等 多 
种 场景 。 特 征 子 系统 保留 了 交流 网 络 与 跟 网 型 装备 


I 


fF 


的 动态 特性 ， 存 在 与 之 对 应 的 低 维 物理 系统 ， 故 由 
此 导出 的 广义 短路 比 计 算 和 系统 强度 量化 方法 具 
有 较 严 说 的 理论 基础 ， 也 有 较 清 晰 的 物理 意义 。 

本 文 聚焦 跟 网 型 装备 主导 的 稳定 模 态 ， 其 特征 
子 系统 为 跟 网 型 装备 的 单机 无 穷 大 系统 ; 其它 模 态 
下 ， 如 构 网 型 装备 控制 特性 不 能 忽略 时 ， 特 征 子 系 
统 还 需要 保留 构 网 型 装备 ， 其 特性 需 深入 研究 。 进 
一 步 ， 特 征 子 系统 可 能 是 电力 系统 的 基本 单元 , 包 
含 但 不 限于 单机 无 穷 大 系统 ， 其 基本 理论 和 方法 后 
续 还 需 进 
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电力 系统 技术 规定 一 


附录 A 引 理 1 的 举例 


以 两 维和 矩阵 为 例 说 明 引 理 1。 设 M 与 让 分 别 
为 : 
m=; ilera, Ae (29) 
HH, &ER, UEB, i=1,2,3,4. 
根据 式 (29), 利用 M 的 两 个 特征 子 空间 可 求解 
应 广义 特征 矩阵 的 近似 ; 
L ~Y" MX, 
0370 1,4] 50% ie 
0.752, +0.24£, 
L, ~Y MX, 
=53781,+| 6070 ie 
0.242, +0.76£, 


AH: Ea 0.57] 81, X, =[-042 -091] 81, 
“=[-092 042]®L, Y¥ =[-0.57 -0.84]®L, 


其 中 ， arene 重 特征 值 -0.37 
55.37 的 右 特征 向 量 , Yi. Yo PA el ee At MAE 


特征 向 量 ， ican 
如 果 Ely €2n €3、 647) 


E 子 空间 。 
分 别 为 0.4、0.3、0.2、0.1， 


sth: Ts) -=A()+B eC) Ls)=B (>) 
1,(3) =B,26(s),T,(3) -Be() 
将 子 矩 阵 表 达 式 代入 并 化 简 ， 可 得 ; 
T,(s)-T,(s)T,"(s)T,(s) 
=Y (8) + By, @6(s) 
附录 C 引 理 2 的 证 明 


(35) 


根据 特征 值 与 特征 子 空间 定义 可 知 ， 
(U7/P)B’ rea @ 严 有 二 个 二 重 特征 值 帮 ， 对 应 有 一 组 
特征 子 空间 只 (XY,) 、R(Y) ， 其 中 , X= x; @ be 
Y=y/ @h, Xr YHA A(O/P)B ea Qh 
关于 AWA. ARES. Alt, Mee 1 式 


(3), FAR(X,) AUR(Y,) RAW P Brea Dh tH n 
个 广义 特征 矩阵 ; 
L =¥}'((U?/P)B,,, @L,)X,=24,8L, (36) 
WE, RAR RIE E m PSG OR th AL, 
M, XA Ye Pet: 
MX, =X,L,,Y"M=LY" 
式 中 : L, =Y pra (S )+4, @¢(s ) 
比较 式 (1)、(37) 并 结合 特征 子 空间 定义 (定义 1) 
mga, R(X). R(Y,) 和 工分 别 为 原 系统 矩阵 的 特 
征 子 空间 及 广义 特征 和 矩阵, 因此 , 令 X=(Xq... Xna) 
Y=(Yq... Yng)s 利用 X、 了 可 将 原 系统 矩阵 转换 为 : 
L=Y"MX =diag(L,) (38) 
BFAA Y FAS EREJE, We det (M1) =0 
5 det(L)=0 Ft, BS det(L,)=0% fr. AL, X 


应 的 特征 子 系统 (上 为 特征 子 系统 矩阵 ) 可 近似 原 
系统 的 稳定 性 。 证 毕 。 
附录 D 定理 1 的 证 明 
TERA: 利用 只 (XX,) . RO, ) 可 求解 等 效 同 构 系 
SFE ME Mi 的 广义 特征 矩阵， 满足 : 
L, =Y} M,X, =G,(s)+A,-¢(s) 


(37) 


(39) 


第 刘 晨 曦 等 ， 基于 特征 子 系统 的 广义 短路 比 导 出 原理 及 计算 方法 13 
AE Dy V AESTAE 7 
结合 式 (39)、(19), 等 效 同 构 系 统 的 广义 特征 矩 9 人 馈 入 系统 电网 参数 (p.u.) 
E v AET $ ` H X15=0.20 X2, 6=0.40 X3,7=0.60 Xi 2=0.40 X23=0.40 
年 是 原 系统 广义 特征 矩阵 的 近似 。 根 据 引 理 1， 利 i J a 2 z 
X13=0.40 X3,4=0.20 
MIAE V yyt = Ap YE Te Jer 7 7 
J R(X, ) 和 R(Y ) TE SUREE EAD. HEHE ELGG ESR 0 pu.) 
Xss,86=0.24 Xs6.86=1.00 X57,36=2.40 X58.86=3.30 X59,86=2.90 
附录 E X60,86=2.00 X61,86=1.70 X62,86=0.03 X63,87=1.10 X64,87=3.10 
EL 新 能 源 及 储 能 变 流 器 控制 参数 ee A E E 
X70,87=0.02  X7,88=2.50 X72,88=2.20 X73,88=1.50  X74,88=1.30 
Table El Control parameters of grid-following converters , i , > 
有 m 储 能 变 流 器 X75,88=0.81 X76.88=1 .30 X77,88=0.81 X78.90=0.19 X79,90=2.20 
内 容 新 能 源 变 流 器 Xgog9=0.19 Xs1,89=2.20 Xg2,91=1.50 Xg3,91=1.30 Xs4,91=1.50 
RR RAE 1500kVA 620V Xss,91=1.30  Xg6s,87=1.00  Xs792=3.10 Xss, 92=4.40 X9077.00 
量 、 电 压 Xoo9%=7.80  X91,92=0.87 Xo2093=6.90 变压器 电抗 : 64.0 
装 名 额定 容量 ENOR 0.68MVA(14,15) 算 例 1 在 gSCR=SCRo 时 的 电网 所 有 参数 是 上 述 电 网 参数 的 1.1 信 
0.72MVA(28,29) 
直流 侧 电 容 、 电 压 0.1089p.u., 1.1kV 0.1089p.u., 1.5kV 表 E4 某 实际 新 能 源 基 地 4 个 算 例 中 的 装备 实际 工 况 
外 环 PI 0.5,40( 直 流 电 压 ) 0.5,20( 功 率 ) Table E4 The operating conditions of devices in four examples 
滤波 电感 、 电 容 0.05、0.0006 0.05、0.0006 of Practical Renewable Energy Plant 
ee 0.5,15(14,15) 算 例 运行 工 况 
电流 内 环 PI 0.2,60 0.3,10(28,29) P 1-13:0.7;16-27:0.8;30-38:0.6;39-42:0.75; 44-53:0.9;14-15: - 
ya m2 ; i 0.45;28-29: -0.48;43,54:1.00; 
馈 入 系统 : 14,10400 yy | 1:13, 16-27, 30-41, 44-48 :1.01;14-15, 28-29:1.03;42,43,54:1.00; 
某 实际 新 能 源 场 站 : 49-53:1.02 
P 1.00 
: 0,10500(1-13); 21,10300 
锁 相 环 PI ms f i 40,24000 ? [v 1.00 
(16-27); 26,10400 (30-38); ie [Se 14-15: 0.45:28-29: 0.48， 其 余 见 算 例 1 
6,10200(39-42,44-48); U 14-15, 28-29:1.01， 其 余 见 算 例 1 
31,10200(49-53) p P 见 算 例 1 
电压 前 馈 时 间 常 数 0.00001 U 见 算 例 1 
表 E2 构 网 型 变 流 器 控制 参数 收 稿 日 期 ; 
Table E3 Control parameters of grid-forming converter ROS 
作者 简介 : 
内 容 数值 JIT : 
装备 额定 容量 1.5MVA 辛 焕 海 (1981)， 男 ’ 教授 ， 通信 作者 ， 博士 生 导 师 ， 主 
滤波 电感 、 滤 波 电容 、 电 阻 0.05、0.0005、0.02 要 研究 方向 为 : 新 能 源 电 力 系统 稳定 分 析 与 控制 。E-mail 
电压 外 环 、 电 流 内 环 PI 参数 0.8,12; 6,10 


构 网 型 变 流 器 惯量 系数 、 阻 尼 系 数 2,20 
电压 前 馈 滤 波 时 间 常 数 0.001 


co 


RE 仿真 电网 参数 


Table E3 Parameters of Simulations 
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Derivation Principle and Calculation Method of Generalized Short-circuit Ratio Based on 


Eigen-Subsystems 


XIN Huanhai!, LIU Chenxi!, HUANG Linbin!, GAO Huisheng!, JU Ping! 
(1. College of Electrical Engineering, Zhejiang University) 


KEY WORDS: generalized short-circuit ratio; device critical short-circuit ratio; eigen-subsystems 


With the large-scale integration of power electronic 
devices, voltage support strength(referred to as “system 
strength”) decreases, thus threatening the system's safe 
and stable operation. In both the homogeneous and heter- 
ogeneous scenario (i.e., power electronic devices inte- 
grated into the system have the same dynamics or weakly 
different dynamics), system strength can be quantified by 
the gSCR and SCRo. However, there is a lack of unified 
derivation principles and methods. To this end, this paper 
unified these based on the concept of eigen-subsystems. 

Specifically, we consider a power system integrated 
with types of devices, where n devices are grid-following 
converters or conventional energy storage converters; m 
devices are grid-forming devices. In the small-signal sta- 
bility analysis, grid-forming devices can be equivalent to 
an ideal voltage source with the series inductor, whose 
equivalent admittance is B1BR (j=1,..., m). Moreover, con- 
sider decoupling the actual operating conditions of the re- 
maining n devices from their dynamic. Thus, the inverse 
matrix of the system closed-loop transfer function is 

M =Y sr (S)+¥ a(S) a) 


where Y (s) = diag(U} /P)B 4 @¢(s) 


Vine (s ) =diag (Fass (s )) 
where Bea represents the Thevenin equivalent admittance 
matrix only containing the n buses and considering the 
equivalent admittance Brsr;; Yrsri(s) is the admittance 
model of the i-th converter, i=1,...,n; P; and U; are power 
output and terminal voltage of the i-th converter. 


An equivalent system matrix M; can be constructed 
in (2), whose eigen-subspace is approximate to its of the 
actual system matrix M since M can be treated as a 
perturbation of M; based on the perturbed theory of eigen- 


subspace. Moreover, the eigen-subspace R(X,) and 


R(Y,) of the 
diag(U?/P;) Brea Q h with respect to its own eigenvalues 


equivalent admittance matrix 
Ais the same as matrix Mi. Therefore, using (X,) and 
R(Y, ) (KE {1,...,n}), a first-order approximation of the 
generalized eigen-matrix can be obtained as (3). Let E; 

be the closed-loop transfer function of the eigen-subsys- 
tem, the eigen-subsystem can approximate the stability of 
the actual system in the corresponding mode. The physical 
meaning of eigen-subsystems is a decoupled subsystem of 
the actual system since it physically retains the dynamic 


characteristics of converter and networks 
M, =1, 8G, (s)+Y,., (s) 


L, ~Y MX, =G,(s)+A,-(s) 
G, (s) =Y” | diag (U? /P)Y nri (s))] X, 
=X ps (ws (5), Pu = Xi Yri 


where X= x:8®h, Vi= yr@h, the column vectors of Xx 
and Y, are the right and left eigenvectors of 
diag(U?/P;)Brea ® D with respect to As; xu and ye are the 
i-th elements of x; and yx, respectively. 


(2) 
(3) 


where 


On this basis, in small-signal stability mode, the sys- 
tem strength depends on the gSCR of the power grid, 
which is defined as the equivalent admittance correspond- 
ing to the eigen-subspace under the 1-th mode of interest: 


gSCR= I"diag(U; /P) BX, 
=1:., (diag(U?/P)B,..) 


where Afm represents solving for the minimum posi- 


(4) 


tive eigenvalue of matrices. 

When the dynamic G,(s) of the equivalent con- 
verter remains unchanged, the smaller the gSCR, the 
weaker the system strength. Based on this monotonic 
characteristic, the minimum gSCR required for the equiv- 
alent converter to reach a critical stable state can be deter- 
mined as (5). SCRo is precisely equal to the critical value 
of the gSCR(CgSCRo). 


SCR, = arg {det(G, (s,)+ gSCR-$(s,))=0} (5) 


gSCR 
According to (6) and (7), it can be determined 
whether the system is stable and whether the system meets 
the stability margin requirements. 


gSCR > CeSCR, = SCR, (6) 
SCR -SC 
es 10% > Po (7) 


The effectiveness of the system strength quantifica- 
tion method is verified in cases. 


